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ABSTRACT 



<n : 

(^) Aims. We have measured the deuterium fractionation and the CO depletion factor (ratio between expected and observed CO abundance) in a 
£■*»•. ■ sample of high-mass protostellar candidates, in order to understand whether the earliest evolutionary stages of high-mass stars have chemical 
00 ' characteristics similar to those of low-mass ones. In fact, it has been found that low-mass starless cores on the verge of star formation have 

, large values both of the column density ratio N(NtD + )/N(N'>Il + ) and of the CO depletion factor. 
Methods. We have observed with the IRAM-30m telescope and the ICMT two rotational lines of N 2 H + and N 2 D + , the (2-1) line of C 17 and 

' DCO + , and the sub-millimeter continuum towards a sample of 10 high-mass protostellar candidates. 

, Results. We have detected N 2 D + emission in 7 of the 10 sources of our sample, and found an average value Af(N?D + )/./V(N 2 H + ) ~ 0.015. This 
O |' value is ~ 3 orders of magnitude larger than the interstellar D/H ratio, indicating the presence of cold and dense gas, in which the physical- 
' chemical conditions are similar to those observed in low-mass pre-stellar cores. Also, the integrated CO depletion factors show that in the 
>| . majority of the sources the expected CO abundances are larger than the observed values, with a median ratio of 3.2. 

Conclusions. In principle, the cold gas that gives origin to the N 2 D + emission can be the remnant of the massive molecular core in which the 
(•--j high-mass (proto-)star was born, not yet heated up by the central object. If so, our results indicate that the chemical properties of the clouds in 
• • ■ which high-mass stars are born are similar to their low-mass counterparts. Alternatively, this cold gas can be located into one (or more) starless 
' core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot decide which is the correct scenario. 

X 

5_l Key words. Stars: formation - Radio lines: ISM - ISM: molecules, continuum 
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1. Introduction 

The initial conditions of the star formation process are still 
poorly understood. In recent years, studies of starless low- 
mass cores have begun to unveil the physical and chem- 
ical features that lead to the formation of low-mass stars 
(Kuiper et al. fl996l Caselli et al. IT9991 Evans et al. I200T1 
Caselli et al. I2002al I2002bl Bergin et al. l2002l Ta falla e t 
al. l2TMll2TK)4ll2006l Harvey et al. l2003l Crapsi et al. 12005V 
On the other hand, the characterisation of the earliest stages of 
the formation of high-mass stars is more difficult than for low- 
mass objects, given their shorter evolutionary timescales, larger 
distances, and strong interaction with their environments. 

In the last years various authors (Molinari et al. 2000, 2002; 
Sridharan et al. 2002; Beuther et al. 120021 Fontani et al. 2005) 
have performed extensive studies aimed at the identification of 
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precursors of ultracompact (UC) Hn regions, i.e. very young 
(< 10 5 yr), massive (M> 8M ) objects which have not yet 
ionised the surrounding medium. In particular, from CS and 
mm continuum observations, it has been noted that in the ear- 
liest stages prior to the onset of massive star formation, the 
radial distribution of the intensity is quite flat, resembling the 
structure of starless cores in more quiescent and less massive 
molecular clouds (Beuther et al. 2002). More evolved objects 
show more centrally-peaked density structures, with power-law 
indices once again resembling those found in low-mass cores 
(i.e. n Hn _oc r' 1 - 6 ; Shu 1977; Motte & Andre 2001). These results 
strongly suggest to apply the investigative techniques that are 
successful in the study of low-mass star forming cores, to the 
high-mass regime. 

It has been found that when the starless core is on the verge 
of dynamical collapse, most of the high density tracers, in- 
cluding CS and the CO isotopologues, are frozen onto dust 
grains. CS observations of low-mass cores clearly show that 
this molecule is avoiding the central high density core nucleus, 
where the mm continuum peaks (Tafalla et al. 2002). The mor- 



2 



Fontani et al.: Massive pre-stellar cores 



phology of the C n O (1-0) and (2-1) line emission of IRAS 
23385+6053, a reliable example of precursor of an UC Hn re- 
gion (Fontani et al. 2004a), suggests that CO is depleted in the 
high-density nucleus traced by the continuum emission. On the 
other hand, species such as NH3, N2H + and N2D + are good 
tracers of the dust continuum emission (Tafalla et al. 2002; 
Caselli et al. 2002b; Crapsi et al. 2004, 2005), suggesting that 
they are not significantly affected by freeze-out (see Bergin et 
al. 2001 and Crapsi et al. 2004 for some evidence of N2H + de- 
pletion toward the center of B68 and L1521F, respectively). 
This is probably due to the fact that the parent species N2, 
unlike CO, maintains a large gas phase abundance at densi- 
ties around 10 6 crrr 3 , despite the recent laboratory findings 
that N2 and CO binding energies and sticking coefficients are 
quite similar (Oberg et al. 2005; Bisschop et al. 2006). The au- 
thors above mentioned also found a relation between the col- 
umn density ratio D frac =N(N2D + )/N(N2H + ) and the core evo- 
lution: Df rac is predicted to increase when the core evolves to- 
wards the onset of star formation, and then it is expected to 
drop when the young stellar object starts to heat up its sur- 
roundings. Summarising, pre-stellar cores or cores associated 
with very young massive stars should show strong emission in 
the N2H + and N2D + lines, with large values of Df rac , and little 
or no emission of molecular species such as CO and CS, if they 
evolve like their low-mass counterparts. 



Although recently Sridharan et al. (2005) and Beltran et 
al. (2006) have detected several candidate high-mass starless 
cores in the neighbourhood of some massive protostar candi- 
date, a sample of well-identified high-mass pre-stellar cores 
is lacking. At present, the objects that are believed to be the 
closest to the earliest stages of the high-mass star formation 
process are those of two samples of high-mass protostellar can- 
didates, selected from the IRAS-PSC by Molinari et al. ( 1996) 
and Sridharan et al. (2002i, and then investigated by various 
authors (Molinari et al. fT998ll2000l Brand et al. l200Tl Fontani 
et al. 2004a ,b; Beuther et al. l200l Williams et al. l2004l Fuller 
et al. 2005 1. These studies have shown that a large fraction of 
these sources are indeed very young massive objects. Also, as- 
suming that high-mass stars form in clusters (Kurtz et al. l2000> . 
we may expect to find other cold and dense cores located close 
to our target sources. For this reason, we have observed with 
the IRAM-30m telescope and SCUBA at JCMT a selected sam- 
ple of candidate precursors of UC Hn regions in N2FU, N2D + , 
C 17 and sub-mm continuum. The sources have been selected 
from the two samples mentioned above on the basis of observa- 
tional features indicative of very little evolution: they are asso- 
ciated with massive (M ~ 1O 2 M ) and compact (diameters less 
than 0. 1 pc) molecular cores, in which the kinetic temperatures 
are around ~ 20 K, and show faint or no emission in the radio 
continuum. Also, they have distances less than 5 kpc. 



Sect. 13 describes the observations and the data reduction, 
while the results are presented in Sect. and discussed in 
Sect. [4] A summary of the main findings is given in Sect. [5] 



2. Observations and data reduction 

2.1. Molecular lines 

All molecular tracers were observed with the IRAM-30m 
telescope. In Table ^ we gi ve the molecular transitions ob- 
served (Col. 1), the line rest frequencies (Col. 2), the telescope 
half -power beam width (HPBW, Col. 3), the channel spacing 
(Col. 4) and the total bandwidth (Col. 5) of the spectrometer 
used. The observations were made in wobbler-switching mode. 
Pointing was checked every hour. The data were calibrated with 
the chopper wheel technique (see Kutner & Ulich 1981 1, with 
a calibration uncertainty of ~ 20%. 

2.1.1. N 2 H + and N 2 D + 

A sample of 19 sources, selected from two sample of proto- 
star candidates as explained in Sect. 1, were observed in July 
2002, obtaining single-point spectra of the N2H + (1-0) and (3- 
2) lines, and N2D + (2-1) and (3-2) lines towards the position 
of the IRAS source. Then, according to the results obtained, in 
a second observing run carried out in July 2003, we mapped 
in on-the-fly mode the 10 brightest sources in both the N2FU 
transitions, in order to have a clear identification of the emis- 
sion peak position. Finally, spectra of the N2D + (2-1) and (3-2) 
lines were obtained with a deep integration towards the N2FU 
(1-0) line peak position. These sources are listed in Table |2] 
while those observed during the first run only are listed in 
Table|3] 

The on-the-fly maps of the N2FU (1-0) and (3-2) transi- 
tions were obtained in a 1' xl' field centered on the position 
listed in Table|2] The two transitions were mapped at the same 
time and registered with two spectrometers: one with low spec- 
tral resolution and large bandwidth and a second with higher 
spectral resolution and narrower bandwidth (see Table^. The 
antenna temperature, 7*!, and the main beam brightness tem- 
perature, 7mb are related as T* A = r M B ?7mb, with 77MB = 
fl e ff/^eff=0.80 and 0.48 for the (1-0) and (3-2) lines, respec- 
tively. The spectra of the N2D + (2-1) and (3-2) lines were ob- 
tained simultaneously, again with two spectrometers with dif- 
ferent spectral resolutions and bandwidths (see Table [Q. The 
values of 77MB are 0.73 and 0.57 for the (2-1) and (3-2) transi- 
tions, respectively. 

The N 2 H + (1-0), (3-2), and N 2 D + (2-1), (3-2) rotational 
transitions have hyperfine structure. To take this into account, 
we fitted the lines using METHOD HFS of the CLASS pro- 
gram, which is part of the GAG software developed at the 
IRAM and the Observatoire de Grenoble. This method assumes 
that all the hyperfine components have the same excitation tem- 
perature and width, and that their separation is fixed to the lab- 
oratory value. The method also provides an estimate of the total 
optical depth of the lines, based on the intensity ratio of the dif- 
ferent hyperfine components. The frequency of the N2FU (1- 
0) line given in Tabled is that of the main component (Dore 
et al. I2004> . while that of the (3-2) transition has been taken 
from Crapsi et al. (2005 1, and it refers to the F\F - A 5 — > 3 4 
hyperfine component, which has a relative intensity of 17.46%. 
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Rest frequencies of the N2D + lines given in Table^refer to the 
main hyperfine component (Dore et al. 2004). 

2.1.2. C 17 Oand DCO + 

On-the-fly maps of the C n O (2-1) and DCO + (2-1) transi- 
tions were obtained on July 2003. These two transitions and 
the N2H + (1-0) and (3-2) lines were mapped at the same time 
and with two spectrometers. The data were reduced using the 
CLASS program. Like the N 2 H + lines, the C n O (2-1) line has 
hyperfine structure. However, the spectra were too noisy to re- 
solve the different components, so that the lines were fitted with 
a single Gaussian. 

2.2. Sub-millimeter continuum 

Continuum images at 850 fim of 5 of the 10 sources mapped in 
N 2 H + were taken on 1998 October with SCUBA at the JCMT 
(Holland et al. I1998i . The standard 64-points jiggle map ob- 
serving mode was used, with a chopper throw of 2 ' in the 
SE direction. Telescope focus and pointing were checked using 
Uranus and the data were calibrated following standard recipes 
as in the SCUBA User Manual (SURF). One of the maps of this 
dataset, that of IRAS 22172+5549, has been already published 
by Fontani et al. J2004b> . 

3. Results 

Parameters and detection statistics of the 10 sources ob- 
served more extensively (see Sect. 12.1.11 are listed in Table|5] 
Column 1 gives the IRAS name, and the equatorial (J2000) 
coordinates of the maps center are listed in Cols. 2 and 3. In 
Cols. 4 and 5 we give the source velocities and kinematic dis- 
tances, respectively. Cols. 6, 7, 8, 9 and 10 give the follow- 
ing information: detection (Y) or non-detection (N) in N 2 H + , 
N 2 D + , C n O, DCO + and millimeter or sub-millimeter contin- 
uum, respectively. In Col. 6 we also give the offset of the posi- 
tion of the N 2 H + (1-0) line emission peak with respect to the 
map center. As explained in Sect. 12. 1.11 the N 2 D + spectra have 
been taken toward this position. For the lines mapped, we have 
considered as detected those sources showing emission above 
the 3<t level in the map field. Main parameters of the 9 sources 
observed only in the first observing run are listed in Table|5J the 
equatorial coordinates in Cols. 2 and 3 represent the observed 
position, while source velocity, Vlsr, and kinematic distance, 
ofkin, are listed in Cols. 4 and 5, respectively. In Cols. 6 - 9 we 
also give the peak temperature of the lines detected, or the rms 
level of the spectrum for those undetected. 

In the following, we will consider only the results obtained 
for the sources listed in Table |2] Since they have been previ- 
ously observed in different tracers by other authors, for com- 
pleteness in Table |4] we give the most important physical pa- 
rameters reported in the literature: the dust temperatures T^, 
obtained from gray-body fits to the observed spectral energy 
distribution of the source, the gas kinetic temperature 7\; n , de- 
rived from NH3 observations (Molinari et al. 1996 Jijina et 
al. [19991 or CH 3 C 2 H observations (Brand et al. l200lll . and H 2 



column- and volume densities computed from continuum ob- 
servations. 

3.1. A/ 2 /-/ + and N 2 D + spectra at the N 2 H + emission 
peak 

In the following, we will discuss the data taken with the spec- 
trometer with the highest spectral resolution. All sources listed 
in Table [2] have been detected both in the N 2 H + (1-0) and 
(3-2) transitions. Five have also been detected in both N 2 D + 
(2-1) and (3-2) lines, and two (IRAS 18511+0146 and IRAS 
22172+5549) in the N 2 D + (2-1) transition only. All spectra 
of the N 2 H + (1-0) and N 2 D + (2-1) lines are shown in the 
Appendix, in Figs. lA-l1 and lA-2l 

In Table|5]we give the N 2 H + (1-0) and (3-2) line parame- 
ters of the spectra taken at the peak position of the (1-0) line 
emission: in Cols. 3 - 7 we list integrated intensity (J Tmb^v), 
peak velocity (Vlsr), FWHM, opacity of the main component 
(tio), and excitation temperature (r ex ) of the N 2 H + (1-0) line, 
respectively. Cols. 9-12 show the same parameters for the 
N 2 H + (3-2) line, for which we do not list T ex because we will 
not use it in the following. The integrated intensities have been 
computed over the velocity range given in Cols. 2 and 8, while 
for the other parameters we have adopted the fitting procedure 
described in Sect. 12.1. 11 This procedure cannot provide r ex for 
optically thin lines. However, one can note that for the sources 
with opacity > 0.4, the values of T ex are comparable to the ki- 
netic temperature T^ n listed in Col. 5 of Table0] within a factor 
~ 2 (with the exception of IRAS 19092+0841, for which the 
difference is a factor 4). Therefore, for sources with optically 
thin lines, we decided to assume r ex =7\j n . 

Most of the hyperfine components of the N 2 H + lines are 
not well resolved. This is due to the fact that the line widths 
(-1-3 km s _1 ) are larger than the separation in velocity of 
most of the different components. As an example, we show in 
Fig. [U the spectra of IRAS 05345+3157 with the position of 
the hyperfine components. Regarding the assumption of equal 
excitation conditions made to treat the hyperfine structure of 
the N 2 H + (1-0) line, Daniel et al. J2006> have recently found 
with theoretical calculations that this assumption is not valid in 
clouds with densities in the range 10 4 — 10 6 crrT 3 and line opti- 
cal depths larger than 20. Since all of our sources have compa- 
rable densities (see Sects. I3~3l and l3.4> but the lines have much 
smaller opacities (see Table|5}, the assumption of equal excita- 
tion conditions is considered to be valid. 

The line parameters of the N 2 D + transitions, derived us- 
ing the same fitting procedure adopted for N 2 H + , are listed in 
Table |6] We give an upper limit for the integrated intensity of 
the lines for the undetected sources, using the observed rms and 
an average value of the line FWHM of 1.5 km s , both for the 
(2-1) and (3-2) transitions. With respect to Table|5] we do not 
give the line opacities and T ex because in most of the spectra 
the uncertainties are comparable to or larger than the values 
obtained. The two sources with well-determined opacity of the 
main component of the N 2 D + (2-1) line are IRAS 20126+4104 
and IRAS 20216+4107, for which we obtain t = 0.13 ± 0.09 
and 0.96+0.03, respectively. For the remaining ones, i.e. IRAS 
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Table 1. Observed transitions 



molecular 


frequency 


HPBW 


Av T 


Bandwidth 1 


transition 


(GHz) 


(") 


(km s-') 


(km s~') 


N 2 H + (1-0) 


93.173772 


26 


0.031/3.218 


113/824 


N 2 H + (3-2) 


279.511863 


9 


0.042/1.073 


75/275 


N 2 D + (2-1) 


154.217137 


16 


0.038/1.944 


68/498 


N 2 D + (3-2) 


231.321966 


10 


0.051/1.296 


91/332 


C n O(2-l) 


224.714385 


11 


0.052/1.334 


93/342 


DCO + (2-1) 


144.077319 


16 


0.041/2.081 


73/533 



the two values refer to the two spectrometers used. 



Table 2. Source list and detection summary 



IRAS 


R.A. (J2000) 


Dec. (J2000) 


^LSR 


"kin 


N 2 H + f 


N 2 D + 


c"o 


DCO+ 


(sub-)mm 


Source 


(h m s) 


(o / //) 


(km s" 1 ) 


(kpc) 












IRAS 05345+3 157 1 


05 37 52.4 


32 00 06 


-18.4 


1.8 


Y (0-5) 


Y 


Y 


Y 


Y 


IRAS 05373+2349 1 


05 40 24.5 


23 50 55 


2.3 


1.2 


Y (-8,0) 


Y 


N 


Y 


Y 


IRAS 18144-1723 1 


18 40 24.5 


-17 22 13 


47.3 


4.3 


Y (-7,0) 


N 


Y(?) 


N 


Y 


IRAS 18511+0146 1 


18 53 38.1 


01 50 27 


56.8 


3.9 


Y (0,-5) 


Y* 


Y 


Y 


Y 


IRAS 18517+0437 2 


18 54 13.8 


04 41 32 


43.9 


2.9 


Y (12,0) 


N 


Y 


N 


Y" 


IRAS 19092+0841 1 


19 11 38.6 


08 46 30 


58.0 


4.5 


Y (8,-4) 


N 


Y 


N 


Y 


IRAS 20126+4104 2 


20 14 26.0 


41 13 32 


-3.8 


1.7 


Y (-5,0) 


Y 


Y 


Y(?) 




IRAS 2021 6+4 107 2 


20 23 23.8 


41 17 40 


-2.0 


1.7 


Y(-5,5) 


Y 


Y 


N 




IRAS 20343+4129 2 


20 36 07.1 


41 40 01 


11.5 


1.4 


Y (-8,4) (13,0) 


Y 


Y 


Y(?) 




IRAS 22172+5549 1 


22 19 09.0 


56 04 59 


-43.8 


2.4 


Y (-5,0) 


Y* 


Y 


N 


Y 



1 = Vlsr and from NH3 observations (Molinari et al. 119961 

2 - Vlsr ar >d d kin from NH 3 observations (Sridharan et al. 12002 > 

" f = between brackets, the offset (in arcsec) of the N 2 H + (1-0) line emission peak with respect to the map center (coordinates in Cols. 2 and 3) 
is given 

* = detected in the N 2 D + (2-1) transition only 
(?) = marginal detection 

" = observed by Beuther et al. ( 2002 1 at 1.2 mm 

* = observed by Williams et al. 1 2004 1 at 850 /im and/or 450 yum 



05345+3157, IRAS 05373+2349, IRAS 18511+0146, IRAS 
20343+4139 and IRAS 22172+5549, we have fitted the lines 
forcing the optical depth to be 0.1. The observed LSR veloci- 
ties and line widths are in good agreement with those derived 
from N2H + , indicating that the two molecular species are trac- 
ing the same material. 

3.2. Distribution of the integrated intensity 

The maps of the integrated intensity of the lines (N2H + (1- 
0) and (3-2), C 17 (2-1) an d DCO + (2-1)) are shown in the 
Appendix, from Fig. IA-31 to IA-121 For the sources observed 
with SCUBA, the emission map of each line has been superim- 
posed on the continuum map at 850yum, while for the others the 
position of the (sub-)millimeter peaks found in literature are in- 
dicated. We have not shown the maps of the sources marginally 
detected in C n O and DCO + (2-1) lines. 

Figures lA-31 - IA- 12l show that in all sources mapped with 
SCUBA the distribution of the 850/mi emission and the N 2 H + 
(1-0) and (3-2) integrated line emission are in good agreement. 
In fact, the emission contours at half of the maximum (FWHM) 



are fairly well overlapping. For almost all of the sources, the an- 
gular separation between the lines emission peaks and that of 
the sub-mm emission (or the 1.2 mm for IRAS 18517+0437) 
is much smaller than the beam size, with the exception of the 
N 2 H + (3-2) line in IRAS 18144-1723 and IRAS 18517+0437, 
for which the peak of the map is offset by ~ 12" from the 
dust peak. This result confirms that the N 2 H + molecule and 
the dust continuum emission trace similar material, as already 
found by several authors which observed the N 2 H + lines in 
both low-mass and high-mass star formation regions (Caselli 
et al. l2()()2al Crapsi et al. l2005l Fuller et al. 12005V 

The distribution of the C 17 (2-1) line integrated emis- 
sion follows that of N 2 H + and the 850//m continuum in 
IRAS 18511+0146, IRAS 18517+0437, IRAS 19092+0841, 
IRAS 20216+4107 and IRAS 20343+4129, while in IRAS 
22172+5549 it has a different distribution and in the other 
sources the signal is faint and irregular. Finally, we note that 
the integrated intensity of the DCO + (1-0) line, clearly de- 
tected only toward IRAS 05345+3157, IRAS 05373+2349 and 
IRAS 18511+0146, presents significant differences from one 
source to another: in IRAS 18511+0146 (Fig. lA^6l the DCO + 
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Table 3. Sources observed during the first observing run only. Cols. 6-9 give the peak of each line or the rms level of the 
spectrum, both in K. 



IRAS 


RA. (J2000) 


Dec. (J2000) 


Vlsr 




N 2 H + (1-0) 


N 2 H + (3-2) 


N 2 D+ (2-1) 


N 2 D + (3-2) 


Source 


(h m s) 


(o / //) 


(km s" 1 ) 


(kpc) 










IRAS 05490+265 8 2 


05 52 12.9 


26 59 33 


0.8 


2.1 


0.4 


<0.9 


< 0.19 


<0.3 


IRAS 18567+0700 1 


18 59 13.6 


07 04 47 


29.4 


2.16 


0.5 


< 0.4 


< 0.13 


< 0.2 


IRAS 20106+3545 1 


20 12 31.2 


35 54 46 


7.8 


1.6 


0.5 


0.45 


< 0.05 


<0.05 


IRAS 20205+3948 2 


20 22 21.9 


39 58 05 


-1.7 


4.5 


<0.03 


<0.1 


< 0.05 


< 0.05 


IRAS 20278+352 1 1 


20 29 46.9 


35 31 39 


-4.5 


5.0 


0.3 


< 0.2 


< 0.06 


<0.08 


IRAS 22134+5834 2 


22 15 09.1 


58 49 09 


-18.3 


2.6 


1.2 


1.6 


< 0.06 


<0.08 


IRAS 23139+5939 2 


23 16 09.3 


59 55 23 


-44.7 


4.8 


0.9 


1.5 


< 0.16 


<0.3 


IRAS 23151+5912 2 


23 17 21.0 


59 28 49 


-54.4 


5.7 


< 0.08 


<0.5 


< 0.13 


< 0.36 


IRAS 23545+6508 2 


23 57 05.2 


65 25 11 


-18.4 


0.8 


< 0.07 


<0.8 


< 0.15 


< 0.42 



1 = Vlsr and d^in from NH3 observations (Molinari et al. H996i 

2 - Vlsr and d^ m from NH3 observations (Sridharan et al. 2002 



Table 4. Main physical parameters of our sources found in literature: dust temperature (T^), H2 column and volume densities 
(A^h, and n^, respectively) and gas kinetic temperature (7\j n ). The references are shown between brackets. 



Source 


Id 


Mj 2 


«H 2 


Tkin 




(K) 


(xlO 22 cm~ 2 ) 


(XlO 5 cm" 3 ) 


(K) 


05345+3157 




8.6(1) 


1.5(1) 


17.0 (2) 


05373+2349 


27(1) 


5.9(1) 


2.4(1) 


21.2 (3) 


18144+3157 


29(1) 


17.7 (1) 


1.8(1) 


23.6 (3) 


18511+0146 


35(1) 


6.2(1) 


0.4(1) 


27.3 (4) 


18517+0437 


38 (5) 


68 (6) 


5.4 (6) 


29 (3) 


19092+0841 


29(1) 


12.4(1) 


1.6(1) 


40.4 (4) 


20126+4104 


62(5) 


52 (6) 


11.4 (6) 


23.0 (5) 


20216+4107 


46(5) 


18(6) 


3.1 (6) 


21.0 (5) 


20343+4129 


44(5) 


22 (6) 


3.9 (6) 


18(5) 


22172+5549 


30(1) 


4.8 (1) 


0.8(1) 


17.5 (3) 



(1) Molinari et al. i2000l 

(2) Jijina et al. fP%9l 

(3) Molinari et al. il996l 

(4) Brand et al. <200T1 

(5) Sridharan et al. <2"0U2l 

(6) Beuther et al. <2002l 



emission is in good agreement with that of the continuum; in 
IRAS 05345+3157 it is more extended (see Fig. lA-3> . while in 
IRAS 05373+2349 the two tracers are almost totally separated 

(Fig. EH}. 

Assuming that the emission in each tracer has a Gaussian 
profile, we have derived the angular diameter (9) of the sources 
deconvolving the observed FWHM in that tracer with the ap- 
propriate Gaussian beam. The results are listed in Table 
We could not compute the diameters of IRAS 19092+0841 
and IRAS 20343+4129 in the N 2 H + (1-0) line and of IRAS 
22172+5549 in the N2H + (3-2) line because these are not re- 
solved. The diameters derived from the N2H + (1-0) line range 
from 18.1" to 41.3", while those from the (3-2) transition 
are between 14.9" and 29.9", indicating that the higher ex- 
citation transition traces a more compact region. The values 
of 9 derived from C n O are typically in between those de- 



duced form the N2H + (1-0) and (3-2) lines, with the excep- 
tion of IRAS 05345+3157 and IRAS 18517+0437. However, 
one has to keep in mind that 9 has been derived under the hy- 
pothesis of Gaussian source, which is only a rough assump- 
tion for the C n O emission maps. For the sources mapped with 
SCUBA by Williams et al. (I2U041 i.e. IRAS 20126+4104, 
IRAS 20216+4107 and IRAS 20343+4129, the authors do not 
give any estimate of the angular diameter. For this reason, in 
Table[7]we list the diameters estimated from 1 .2 mm continuum 
maps, derived calculating the geometric mean of the major and 
minor axes of the sources, obtained by Beuther et al. (2002i 
from two-dimensional Gaussian fits (see their Table 2). 

Incidentally, we note that the diameters derived from the 
850/im continuum emission are on average smaller than those 
derived from the 1.2 mm continuum by Beuther et al. (|5002i. 
Since the angular resolutions are comparable, we believe that 



Table 5. N 2 H + line parameters 











N 2 H + (1- 


-0) 










N 2 H + (3-2) 






u 
u 


source 


vel. range 


J r M Bdv 


VLSR 


FWHM 


Tio 


T 

* ex 


vel. range 


J r M isdv 


VLSR 


FWHM 


T32 


o 
o 

M 




(kms" 1 ) 


(Kkms- 1 ) 


(kms- 1 ) 


(kms" 1 ) 




(K) 


(kms- 1 ) 


(Kkms- 1 ) 


(km s" 1 ) 


(kms" 1 ) 




05345+3157 


-28,-9 


14.66 


-17.97 


1.59 


0.5±0.2 


22.0 


-23,-13 


16.33 


-18.03 


1.75 


4.3+0.5 




05373+2349 


-8,11 


20.86 


2.31 


1.57 


0.8+0.3 


20.14 


-2,7 


26.13 


2.26 


1.47 


3.7+0.2 


A 


18144-1723 


36,60 


38.44 


48.01 


3.47 


0.49±0.03 


26.23 


42,53 


31.41 


48.42 


3.01 


7+1 




18511+0146 


46,68 


26.76 


57.12 


2.86 


0.1 




52,62 


22.11 


56.91 


3.33 


1+1 


> 


18517+0437 


33,54 


26.01 


43.84 


2.72 


0.1 




39,48 


28.38 


43.71 


3.21 


0.7+0.3 


s 

■is 


19092+0841 


45,70 


13.20 


57.64 


3.76 


0.5+0.1 


10 


50,63 


13.13 


57.07 


5.84 


0.1 


20126+4104 


-14,7 


36.15 


-3.87 


2.01 


0.1 




-8,1 


39.92 


-3.36 


1.57 


7.82+0.04 




20216+4107 


-12,7 


9.39 


-1.66 


1.30 


1.2+0.3 


9.4 


-6,2 


14.41 


-1.89 


0.99 


8.6+0.5 


20343+4139 a 


0,22 


14.22 


10.87 


1.43 


1.79±0.07 


9.2 


7,14 


14.45 


10.81 


1.21 


0.1 


1 


20343+4139" 


0,22 


18.39 


11.48 


2.44 


1.37+0.07 


9.3 


7,14 


21.51 


11.75 


3.69 


0.1 




22172+5549 


-54,-34 


6.97 


-43.59 


2.04 


0.1 




-49,-39 


14.65 


-43.56 


2.63 


0.1 



a = peak at (-8" ,4") 
b = peakat(13",0") 
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Fig. 1. Spectra of IRAS 05345+3157 in N 2 H + (1-0) and (3-2) lines (left panel) and N 2 D + (2-1) and (3-2) lines (right panel). 
The lines under the N 2 H + spectra indicate the position of the hyperfine components. 

Table 6. N 2 D + line parameters 







N 2 D + (2-1) 






N 2 D + (3-2) 




source 


vel. range 

(km s" 1 ) 


J r M Bdv 

(Kkms" 1 ) 


VLSR 

(km s" 1 ) 


FWHM 

(km s- 1 ) 


vel. range 
(km s" 1 ) 


J 7/MBdv 

(Kkms" 1 ) 


VLSR 

(km s" 1 ) 


FWHM 

(km s-') 


05345+3157 


-19,-16 


0.33 


-17.68 


1.01 


-20,-15 


0.44 


-17.70 


1.57 


05373+2349 


0,4 


0.30 


2.43 


1.65 


0,7 


0.43 


2.38 


1.78 


18144-1723 




<0.45 








<0.68 






18511+0146 


55,60 


0.23 


56.90 


0.99 




<0.54 






18517+0437 




<0.60 








<0.42 






19092+0841 




<0.60 








<0.68 






20126+4104 


-6,0 


0.49 


-4.14 


1.52 


-6,-1 


0.31 


-3.86 


1.05 


20216+4107 


-4,1 


0.31 


-1.57 


1.06 


-4,6 


0.25 


-11.85 


0.72 


20343+4139" 


0,15 


0.48 


10.57 


0.80 


9,14 


0.79 


10.66 


1.21 


20343+4139" 


8,15 


0.60 


11.23 


2.93 


9,15 


0.84 


11.62 


2.17 


22172+5549 


-45,-38 


0.23 


-43.26 


2.2 




<0.14 







: peak at (-8" ,4") 
:peakat(13",0") 



this can be due to the different observation technique. In fact, 
as pointed out in Sect. 12.21 the 850 /im observations have been 
carried out in 'jiggle' map observing mode, which provides 
maps with lower sensitivity with respect to the 1 .2 mm maps 
obtained by Beuther et al. J2002I in the dual-beam on-the-fly 



observing mode, and therefore can be less sensitive to extended 
structures. 
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Table 7. Deconvolved source angular diameters (in arcseconds) derived from each tracer assuming a Gaussian source. 



source 


N 2 H + (l-0) 


N 2 H + (3-2) 


c"o 


(sub-)mm continuum 


DCO + 


05345+3157 


31.8 


20.8 


7.9 


21.4 


34.4 


05373+2349 


34.9 


21.4 




14 


29.6 


18144-1723 


23.8 


22.2 




14.8 




18511+0146 


41.3 


29.9 


41.3 


29.3 


23.6 


18517+0437 


25.2 


18.3 


28.7 


29.0* 




19092+0841 




15.1 


17.6 


11.7 




20126+4104 


32.0 


22.7 


24.7 


18.0* 




20216+4107 


23.6 


14.9 


14.1 


22.7* 




20343+4129 




21.8 


32.5 


26.0* 




22172+5549 


18.1 




24.0 


15.7 





* = from 1.2 mm continuum (Beuther et al. 2002 1. They are the geometric mean of the major and minor axes obtained from two-dimensional 
Gaussian fits to the observed emission. 

IRAS 20126+4104, IRAS 20216+4107 and IRAS 20343+4129 have also been observed at 850 //m by Williams et al. <2004l but they do not 
estimate the diameters. 



3.3. Physical properties from dust emission 

The physical parameters derived form the dust emission are 
presented in Table|8] in Cols. 2 and 3 we give the angular and 
linear diameters (6 cont and D, respectively); in Cols. 4, 5, 6 in- 
tegrated flux density {F v ), gas+dust mass (M cont ) and visual 
extinction (A v ) are listed; Cols. 7 and 8 give the H2 column 
densities (Nu 2 ) derived in two different ways which will be ex- 
plained in the following, and in cols. 9 and 10 the correspond- 
ing H2 volume densities are given. 

The linear diameters have been computed using the kine- 
matic distances listed in Table [2 and are between ~0.08 and 
~0.5 pc, typical of clumps hosting high-mass forming stars 
(Kurtz et al. 2000; Fontani et al.[2005 ). The flux densities, F v , 
are obtained integrating the SCUBA maps by the 3<x level in the 
maps. For IRAS 20126+4104, IRAS 20216+4107 and IRAS 
20343+4129 we give the values listed in Williams et al. (|2004i. 
For IRAS 18517+0437, observed neither in this work nor by 
Williams et al. (2004) at 850 yum, we list the value obtained at 
1 .2 mm by Beuther et al. J2002> . 

From the continuum flux density, assuming constant gas- 
to-dust ratio, optically thin and isothermal conditions, the total 
gas+dust mass is given by: 



which is a plausible 'average' value for high-mass star forma- 



F v d 2 



(1) 



In Eq. Q, d is the distance, k v is the dust opacity coefficient, 
derived according to k Y = k v<j (v/vo'f (where vo = 230 GHz and 
k Vo = 0.005 cm 2 g _1 , which implies a gas-to-dust ratio of 100, 
Kramer et al. H998i . and B V (T) is the Planck function calculated 
at the dust temperature T. For this latter, we decided to use the 
values listed in Col. 5 of Table@](i.e. the gas temperature de- 
rived from NH3 or CH3C2H) instead of the dust temperature 
given in Col. 2 because this latter has been derived from the 
spectral energy distribution of the sources, which includes the 
IRAS measurements, not sensitive to temperatures lower than 
~ 30 K. Therefore, this value can overestimate the dust tem- 
perature throughout the envelope. We have assumed B = 1.5, 



tion regions (see e.g. Molinari et al. 2000, Hatchell et al. 2000 
Williams et al. 2004). 

The visual extinction, A v , has been derived from the con- 



(2) 



tinuum flux density using Eq. (2) of Kramer et al. ( 2003 1: 
1.086 F &50lim 



B&50fim(T)Q. s 



where F^q^ is the integrated flux density at 850/im, and 
Bm)ijm(T) is the Planck function calculated at the temperature 
T, for which we have taken the gas kinetic temperature (Col. 4 
of Table 0. O s is the source solid angle and k' = k^o/kn is 
the ratio of absorption and extinction coefficients. The values 
of Q s have been computed from the angular diameters given 
in Table and we have assumed k' = 1.7 x 10~ 5 , which is 
the average value found by Kramer et al. (2003) for B — 1.5. 
For IRAS 18517+0437, for which F S50lim is not available, we 
scaled the flux at 1 .2 mm to 850 /mi assuming once again a dust 
opacity index B — 1.5 (i.e. a spectral index a = 3.5). 

Finally, we have determined the molecular hydrogen total 
column density following two approaches: (1) from M cont and 
D assuming a spherical source; (2) from A v using the relation 
given in Frerking et al. ( 1982): 



N Hl = 0.94 x 10 21 x A v (mag) [cm" 2 ] 



(3) 



We point out that both approaches use F v but the latter has 
the advantage of not assuming any gas-to-dust ratio, which is 
typically affected by large uncertainties. 

We find masses distributed between ~ 30 to ~ 900 M , 
visual extinctions from ~ 90 to ~ 10 3 magnitudes, average col- 
umn densities between ~ 10 23 and ~ 10 24 ctrT 2 and average 
volume densities of ~ 10 5 — 10 6 cm -3 . Mass, visual extinc- 
tion and gas column density were previously estimated from 
the (sub-)mm continuum emission in six sources of our sample 
by Molinari et al. ( 2000), and in the remaining ones by Beuther 
et al. (2002). Even though both Molinari et al. and Beuther et 
al. used marginally different assumptions in computing the dust 
opacity, and a different approach to determine and A v , their 
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estimates are in good agreement with ours. Table [8] shows that 
the column- and volume density estimates derived from A v is 
~ 3 times smaller than those obtained from M cont . We believe 
that this systematic difference is due to the different hypothe- 
ses made on the dust absorption, since all the other parameters 
used in Eqs. (Q and|2} are the same. However, given the large 
uncertainties associated with these measurements, we conclude 
that the two estimates of Nn^ are in good agreement, indicating 
that the assumed gas-to-dust ratio of 100 is a reasonable value 
for our sources. 

3.4. N 2 H + and N 2 D + column densitites and deuterium 
fractionation 

The deuterium fractionation, Df rac , can be derived by deter- 
mining the ratio of the column density of a hydrogen-bearing 
molecule and its deuterated isotopologue. It is important to 
point out that in cold and dense clouds the molecular species 
used to compute Df rac can be affected by depletion. For ex- 
ample, HCO + and H2CO and their deuterated counterparts are 
highly depleted in the high-density nucleus of low-mass star- 
less cores (e.g. Carey et al. 119981 . Therefore, the deuterium 
fractionation derived from these species is representative only 
of an outer, lower-density shell. On the other hand, it has been 
well established that N2H + and N2D + do not freeze-out even in 
the densest portions of low-mass molecular clouds (e.g. Caselli 
et al. 2002a 2002b , Crapsi et al. 2005 1. Therefore these species 
should give an estimate of the deuterium fractionation repre- 
sentative also of the densest regions of our sources. 

We have derived the N2FF and N2D + column densities, 
jV(N 2 H + ) and jV(N 2 D + ), following the method outlined in the 
Appendix of Caselli et al. (2002b). which assumes a constant 
excitation temperature, T ex . The values of T ex for the N2H + (1- 
0) lines have been derived from the hyperfine fitting procedure 
described in Sect. 12.1.11 For the optically thin lines, for which 
the fitting procedure cannot provide T ex , we have assumed r ex = 
7\i n (Col. 5 of Table|4} as already pointed out in Sect. 13. II For 
the N2D + (2-1) line, for which we do not have good estimates 
of the optical depth, we have assumed the excitation tempera- 
ture of the N 2 H + (1-0) fine. 

The results obtained are listed in Table|5J the N2H + column 
densities are of the order of a few ~ 10 13 cm 2 , while the N2D + 
column densities are a few ~ 10 11 cirT 2 . The values of Df rac 
are between 0.004 and 0.02, with an average value of ~ 0.015. 
This value is ~ 3 orders of magnitude larger than the cosmic 
'average' value of ~ 10~ 5 (Oliveira et al., 2003), and close to 
that found by Crapsi et al. (2005) in their sample of low-mass 
starless cores, indicating that also the physical conditions of the 
gas responsible for the N2D + emission are similar. 

We have also computed the N2H + abundance relative to H2 
by dividing A^(N2H + ) by the molecular hydrogen column den- 
sities listed in Col. 8 of Table[8](i.e. the estimates derived from 
A v ). The average value is ~ 2.05 x 10~ 10 , which is in excellent 
agreement with the 2 x 10~ 10 found by Caselli et al. ( 2002cl in 
a sample of 18 low-mass starless cores. 

Finally, the N2H + and N2D + column densities have been 
computed using the LVG program described in the Appendix 



C of Crapsi et al. ( 2005 1, with which we have also estimated 
the «h, volume densities. The results are given in Table fTol 
We note that both the N2H + and N2D + column densities are 
in good agreement with the estimates given in Table [5] It is 
interesting to compare also the H2 volume density estimates 
given in Tables I8landll0l The values obtained in LVG approx- 
imation are on average ~ 2 times larger than the estimates 
given in Col. 9 of Table [8] and ~ 7 times larger than those 
listed in Col. 10 of Table [S] However, it must be noted that 
the sources with the higher discrepancy are IRAS 1851 1+0146 
and IRAS 18517+0437, for which the density estimates have 
big uncertainties: IRAS 1851 1+0146 is the most distant of our 
sources, while for IRAS 18517+0437 we have derived the flux 
at 850 jt/m from an extrapolation of the 1 .2 mm flux. If we ne- 
glect these two sources, the average ratio between the density 
estimates from Col. 10 of Table [8] and those in Table [Tolbe- 
comes ~ 4. Taking into account that the densities derived in the 
LVG approximation are affected by large uncertainties (even 
1 order of magnitude), we conclude that for the sources with 
more accurate measurements the different density estimates are 
consistent within the uncertainties, confirming that the N2H + 
molecule is a good H2 density tracer. 

Crapsi et al. ( 2005 1 have found that the H2 volume densities 
derived from the LVG code are systematically lower than those 
obtained from dust, suggesting a possible partial depletion of 
N2H + in the inner nucleus of their starless cores. However, 
this nucleus has a diameter of ~2500 A.U. (see e.g. Caselli et 
al. 2005), whose contribution is negligile at a distance of some 
kiloparsecs with our angular resolution. Also, we have to con- 
sider that the H2 volume densities derived from dust assume 
an isothermal cloud, which could not be a good assumption for 
our sources. 

3.5. CO depletion factor 

From the C 17 maps we give an estimate of the integrated CO 
depletion factor, /b, defined as the ratio between the 'expected' 
abundance of CO relative to H2, X~ , and the 'observed' value 



where X° Q is the ratio between the observed CO column den- 
sity and the observed H2 column density. 

The C 17 column density has been derived from the inte- 
grated line intensity at the dust peak, assuming optically thin 
conditions. This assumption is based on the results of several 
surveys of C n O performed in high-mass star formation regions 
with comparable mass and density (e.g. Hofner et al. [2000, 
Fontani et al. 2005). Under this assumption, and using the 
Raleigh- Jeans approximation (valid for our frequencies), one 
can demonstrate that the column density of the upper level, /, 
is related to the integrated line intensity 7(C 17 0) according to: 

Nj = 1.209 10 3 27+1 k , 7(C 17 0) (5) 
/ r\ v [i L v 

where k is the Boltzmann constant, v the line rest frequency, 
T] v the beam filling factor and jj. the molecule's dipole moment 
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Table 8. Source properties derived from (sub-)mm continuum emission: deconvolved angular diameter (6 cont ), integrated 850 
yum continuum flux (F v , for IRAS 18517+0437 we give the integrated 1.2 mm continuum flux), linear size (D), gas+dust mass 
(M cont ), visual extinction (A v ) and H2 column and volume densities (A^h, and «h 2 X derived from M cont and A v . 



source 


$cont 

(arcsec) 


D 

(pc) 


F v 

(Jy) 


M omt 
(M ) 


A v 
(mag) 


N ] ' 

H 2 

(xl0 23 crrT 2 ) 


H 2 

(xl0 23 crrT 2 ) 


«H 2 

(xl0 5 crrT 3 ) 


(xl0 5 crrr 3 ) 


05345+3157 


21.4 


0.19 


4.42 


179 


220 


6.26 


2.07 


10.9 


3.6 


05373+2349 


14 


0.08 


2.87 


35 


239 


6.81 


2.25 


27.8 


9.3 


18144-1723 


14.8 


0.31 


7.76 


1118 


496 


14.1 


4.67 


14.8 


4.9 


18511+0146 


29.3 


0.55 


10.23 


959 


137 


3.89 


1.29 


2.3 


0.8 


18517+0437 


29 


0.41 


7.2 1 


336 


87 


2.47 


0.8 


2.0 


0.6 


19092+0841 


11.7 


0.25 


5.57 


426 


283 


8.05 


2.7 


10.3 


3.4 


20126+4104 


18 


0.15 


21.9 2 


504 


982 


28.0 


9.2 


61.2 


20.2 


20216+4107 


22.7 


0.19 


6.1 2 


160 


196 


5.58 


1.84 


9.7 


3.2 


20343+4129 


26 


0.18 


192 


426 


586 


16.7 


5.51 


30.7 


10.3 


22172+5549 


15.7 


0.18 


3.47 


239 


307 


8.73 


2.88 


15.5 


5.2 



1 = from Beuther et al. 1 2002 1, derived from the 1.2 mm continuum 

2 = from Williams et al. 120041 
'' = from M collt 

t+ = from A v 

Table 9. N 2 H + and N 2 D + total column densities (jV(N 2 H + ) and jV(N 2 D + )), deuterium fractionation (Df rac ), and N 2 H + chemical 
abundance relative to H 2 (Xn,h+)- ^n,h+ has been calculated from the molecular hydrogen column density given in Col. 8 of 
TableEl 



source 


iV(N 2 H + ) 


iV(N,D + ) 




^N 2 H+ 




(xl0 13 cnr 2 ) 


(xlO u cirr 2 ) 




(xio- 10 ) 


05345+3157 


3.14+0.07 


3.2+0.9 


0.010+0.007 


1.5 


05373+2349 


4.34+0.06 


2.8+0.8 


0.006+0.001 


1.9 


18144+3157 


9.75+0.08 


<4.8 


<0.004 


2.1 


18511+0146 


6.47+0.01 


2.5+0.4 


0.004+0.001 


5.0 


18517+0437 


6.48+0.02 


<6.7 


<0.01 


8.1 


19092+0841 


1.9+0.5 


<5.3 


<0.028 


0.7 


20126+4104 


7.62+0.01 


4.9+0.7 


0.006+0.001 


0.8 


20216+4107 


1.4+0.2 


2.8+0.5 


0.020+0.007 


0.8 


20343+4129" 


2.2+0.1 


4.4+0.7 


0.020+0.005 


0.4 


20343+4129" 


2.9+0.1 


5.4+0.8 


0.02+0.01 


0.5 


22172+5549 


1.226+0.006 


2.1+0.4 


0.017+0.004 


0.4 



a = peak at (-8" ,4") 
b = peak at (13" ,0") 



ues obtained from A v . We decided to use these latters rather 
than those derived from M cont because they do not assume any 
gas-to-dust ratio. 

The C 17 'expected' abundance has been obtained for 
each source taking into account the variation of Carbon and 
Oxygen abundances with the distance from the Galactic Center. 
Assuming the standard value of 9.5xl0~ 5 for the abundance 
of the main CO isotopologue in the neighbourhood of the so- 
lar system (Frerking et al. I1982L we have computed the ex- 
pected CO abundance at the Galactocentric distance {Dec) of 
each source according to the relationship: 

X£ = 9.5 10" 5 exp(1.105 - (0.13 D GC (kpc))) , (7) 



(0.1 1 Debye for CO). To compute r\ v , we have used the angular 
diameters listed in Table derived from C n O for all sources 
but IRAS 05345+3157 and IRAS 18144+3157, for which we 
have taken the diameter of the sub-mm continuum. The total 
column density has been obtained from: 

N C "o = — eC^expf-^-) (6) 

where gj and Ej are statistical weight (=27 + 1) and energy 
of the upper level, and Q(T ex ) is the partition function at the 
temperature T ex . We have used as excitation temperature the 
kinetic temperatures listed in Col. 4 of Table @] Then, we have 
computed the observed C n O abundances by dividing A^c"o by 
the H 2 column densities listed in Col. 8 of Table|8] i.e. the val- 
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Table 10. H2 volume densities and N2H + and N2D + column densities, derived in the LVG approximation. 
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N 2 H + 




N 2 D + 


Source 




iV(N 2 H + ) 




A?(N 2 D + ) 




(XlO cm 


-3\ /wiAl3 2\ 

) (XlO cm ) 


(XlO cm 


) (XlO cm ) 


05345+3157 


2.4 


3.2 


8.0 


8.6 


05373+2349 


2.8 


5.0 


9.0 


6.7 


18144+3157 


1.6 


11 






18511+0146 


1.0 


3.9 






18517+0437 


1.7 


3.5 






19092+0841 


1.3 


7.2 






20126+4104 


2.6 


8.6 


6.4 


8.0 


20216+4107 


1.5 


1.1 


1.0 


6.0 


20343+4129" 


0.4 


4.0 


3.0 


2.9 


20343+4129" 


1.0 


9.6 


6.0 


30 


22172+5549 


4.0 


1.4 







a = peak at (-8" ,4") 
b = peakat(13",0") 



which has been derived according to the abundance gradients 
in the Galactic Disk for 12 C/H and 16 0/H listed in Table 1 of 
Wilson & Matteucci ( 1992), and assuming that the Sun has a 
distance of 8.5 kpc from the Galactic Center. Then, follow- 
ing Wilson & Rood ( 1994 1, we have assumed that the oxygen 
isotope ratio 16 0/ 18 depends on Dec according to the rela- 
tionship 16 0/ 18 0=58.8 x £>cc(kpc) + 37.1; finally, taking the 
standard ratio 18 0/ 17 0=3.52 (Frerking et al. 1982 1, we have 
computed the expected abundance of the C n O, X^ llQ , for each 
source according to: 

X E 

X E — (%\ 

c "° 3.52 (58.8 D GC + 37.1) 

The results are listed in Table ITU Cols. 2 and 3 give the 
source Galactocentric distance and the corresponding X E , 
respectively; Cols. 4, 5, 6, 7 and 8 list the C 17 (2-1) in- 
tegrated line intensity (7(C n O)), the beam filling factor (tj v ), 
the C n O column density (N(C 17 0)), the observed C n O abun- 
dance (X£ 17Q ) and the CO depletion factor (f D = X E ll0 /Xg 1Q ). 

Col. 7 of Table ^] shows that in 7 out of 9 sources 
in which we have detected C n O, fu is smaller than 10, 
while in the remaining two sources, IRAS 20126+4104 and 
IRAS 20343+4129, f D is much higher than 10. For IRAS 
1 85 1 1 +0146 we obtain an unusual value of 0.4. We believe that 
this is due to the fact that this source is the most further away 
from the Galactic Center (Dec ~ 12 kpc), and the expected 
CO abundance at such a distance could be different from that 
calculated from Eq. 0. These results indicate that in the ma- 
jority of our sources the observed abundance is comparable to 
or marginally lower than the expected value. However, the dis- 
cussion of these results requires three main comments: first, the 
values of the 'canonical' CO abundance measured by other au- 
thors in different objects varies by a factor of 2 (see e.g. Lacy et 
al l 19941 Alves et al. U999l l. Second, the integrated CO deple- 
tion factor is an average value along the line of sight. Therefore, 
given the large distance to our targets and the poor angular res- 
olution of our data, the effect of not-depleted gas associated 



with the more external gaseous envelope of the source can sig- 
nificantly affect the measured /b, which hence is to be taken 
as a lower limit. Finally, the angular resolution of our observa- 
tions also allow us to derive only average values of fu over the 
sources, which may have complex structure. 

4. Discussion 

In the following we discuss the parameters presented in the pre- 
vious sections, and compare them with the findings from the 
literature. 

4.1. Deuterium fractionation and CO depletion 

Theoretical models predict that the deuterium fractionation is 
correlated to the amount of CO depletion (Caselli et al. 2002b ; 
Aikawa et al. 2005). Observations of low-mass objects have 
partially confirmed this theoretical prediction. Bacmann et 
al. (|2003 ) have found a good correlation between f D and Df rac 
(derived from D2CO/H2CO) in 5 pre-stellar low-mass cores. 
More recently, Crapsi et al. (2005) have obtained Df rac using 
N2DVN2I-F in a sample of 31 low-mass starless cores, and 
they have shown that cores with higher Df rac have also higher 
/d- 

In order to check if such a correlation exists also in high- 
mass objects, we have added to Fig. 5 of Crapsi et al. (2003J, in 
which Df rac is plotted against the CO depletion factor, the val- 
ues obtained from our high-mass objects. The result is shown in 
Fig. |2 we can note that all the objects of our sample have Df rac 
smaller than those of the Crapsi et al. sample, and nearly half of 
them have also smaller fu- The plot also indicates that the cor- 
relation found by Crapsi et al. ( 2005 1 and Bacmann et al. ( 2003 1 
for low-mass sources cannot be extended to the sources of our 
sample. In fact, Df rac does not show any dependence from fy. 

Neverthless, when discussing these results we have to keep 
in mind several caveats. First, the sources of Crapsi et al. ( 2005 1 
are embedded in low-mass molecular clouds ~ 100 - 200 pc 
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Table 11. Parameters used to determine the CO depletion factor: source Galactocentric distance (Dqc), 'expected' C n O abun- 
dance (X^ 17Q ), integrated intensity of C n O at the dust peak position (7(C 17 0))), beam filling factor (r] v ), C n O total column 
density (N(C l1 0)), observed C n O abundance (X°„ ) and CO depletion factor (/ D ). 



source 


D GC 
(kpc) 


K c"o 
(xlO- 8 ) 


/(CO)) 
(Kkms- 1 ) 


1v 


N(C n O) 
(xl0 15 crrr 2 ) 


A c"o 
(xlO" 8 ) 


/d 


05345+3157 


6.38 


8.62 


4.10 


0.34 


6.04 


2.9 


3.0 


05373+2349 


6.96 


7.38 


<2.38 


0.62 


<1.99 


<0.9 


>8.2 


18144+3157 


8.85 


4.62 


6.78 


0.64 


5.61 


1.2 


3.8 


18511+0146 


11.5 


2.56 


12.84 


0.93 


7.73 


6.0 


0.4 


18517+0437 


7.85 


6.0 


8.11 


0.88 


7.39 


9.2 


0.7 


19092+0841 


5.64 


10.6 


9.46 


0.72 


9.03 


3.3 


3.2 


20126+4104 


4.67 


14.3 


5.75 


0.83 


3.65 


0.4 


35.8 


20216+4107 


7.47 


6.47 


7.13 


0.62 


5.93 


3.2 


2.0 


20343+4129 


6.88 


7.54 


3.86 


0.91 


2.17 


0.4 


18.8 


22172+5549 


9.42 


4.05 


2.39 


0.83 


1.45 


0.5 


8.1 



away from the Sun, therefore their estimates of the deuterium 
fractionation and the CO depletion factor, obtained with the 
same angular resolution, are less affected than ours by the con- 
tribution of the non-deuterated and non-depleted gas along the 
line of sight associated with the molecular envelope of the 
source. Additionally, we are comparing a sample of well known 
starless cores, most of which are in the pre-stellar phase, with a 
sample of high-mass protostellar candidates, in which the heat- 
ing produced by the forming protostar may affect both deutera- 
tion and CO depletion. Observations of deuterated molecules 
in low-mass protostars have shown that the deuterium frac- 
tionation is similar to the values found in pre-stellar cores, as 
shown for example by Loinard et al. (2002 1, Hatchell (2003 1, 
and Parise et al. (I2006> . who have observed formaldehyde, 
ammonia and methanol towards several low-mass protostars. 
However, these values are relative to the gaseous envelope in 
which the low-mass protostars are born, which is thought to 
maintain the initial conditions of the star formation process, 
while the heating produced by a forming high-mass protostar is 
expected to dramatically push down the deuteration and the CO 
depletion of the molecular environment (see e.g. Turner 1990 1. 
Finally, as already mentioned in Sect. 13.51 the angular reso- 
lution of the maps allow us to derive only average values of 
Df rac and fo over the sources, which in principle have com- 
plex morphology and may host objects in different evolution- 
ary stages (see e.g Kurtz et al. 2000 1. The origin of such cold 
gas in our sources is thus unclear. We further discuss this point 
in Sect. 1431 

4.2. Effect of temperature on deuterium fractionation 
and CO depletion 

The deuterium fractionation and the CO depletion factor in a 
molecular core are related to its physical properties. In particu- 
lar, in sources with comparable gas density, Df lac and f D should 
be larger in colder clouds. In this section, we will describe a 
simple chemical model which has been used to interpret the 
present observational results on the depletion factor, the deu- 
terium fractionation and the gas temperature. Unlike low mass 



cores, the structure of massive cores cannot be derived in de- 
tail with current observations. The physical properties listed in 
Table 0]and[S]are average values within telescope beams whose 
angular sizes are comparable to source sizes (compare HPBW 
in Table 1 with 6> cont in Table[8}. Therefore, any attempt of mod- 
eling massive cores as smooth objects with densities and tem- 
perature gradients similar to those found in low-mass cores is 
subject to large uncertainties, also considering that massive star 
forming regions are probably clumpy, so that dense and cold 
material may be confined in small regions with filling factors 
significantly smaller than unity. 

For these reasons, we decided to use a chemical model sim- 
pler than the one used for low mass cores (see e.g. Crapsi et al. 
2005). In particular, the detailed physical structure is neglected 
and the clouds are treated as homogeneous objects with den- 
sity and temperature equal to the average values from column 
10 of Tableland column 5 of Table|4] respectively. As already 
pointed out in Sect. 13.31 we used the gas temperature instead 
of the dust temperature to compare with our model predictions, 
given that the latter has been determined using also the IRAS 
measurements, which are not sensitive to temperatures lower 
than 30 K (unlike the data used to determine r gas ), and it is 
thus expected to overestimate the dust temperature throughout 
the cloud. If clumps denser and colder than the average values 
are indeed present, our calculation is expected to underestimate 
the observed CO depletion factor and the deuterium fractiona- 
tion in our sources, unless the filling factor is significantly less 
than unity. 

Other simplifications include the determination of the elec- 
tron fraction, which is now simply assumed to be given by (fol- 
lowing McKee 1989): 

x(e) = 1.3 x 10~ 5 xn(H 2 r°- 5 (9) 

Dust grains (all assumed to be negatively charged l ) are 
distributed in size according to Mathis, Rampl & Nordsieck 

1 Although this is probably true within a factor of 2 (e.g. Flower & 
Pineau des Forets 2003), this assumption does not affect our conclu- 
sions. 
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Fig. 2. Deuterium fractionation (Df rac ) versus integrated CO depletion factor (/b) for our sources (filled circles) and the low-mass 
pre-stellar cores of Crapsi et al. (|2005 1 (open squares). The arrows indicate the upper limits on Df mc (see Table |9]l or the lower 
limit on /d (see Table ll 1> . 



(1977; hereafter MRN) (but with the minimum size value en- 
larged by a factor of 10, as in the case of low-mass cores; if 
we assume the MRN value, conclusions do not appreciately 
change in the range of physical parameters considered here). 
We also included in the model all the multiply deuterated forms 
of and their destruction via dissociative recombination with 
electrons, recombination onto negatively charged grains, and 
through reactions with CO, O and N2 has been taken into ac- 
count. The freeze-out of CO, O and N2 is balanced by thermal 
(see e.g. Hasegawa et al. 1992) and non-thermal (cosmic-ray 
impulsive heating; Hasegawa & Herbst 1993) desorption. In 
analogy with similar work done on low-mass cores, the bind- 
ing energies used are those for CO and N2 onto H2O (see Oberg 
et al. 2005), and the oxygen binding energy is fixed at 800 K 
(Tielens & Allamandola 1987 1. The cosmic-ray ionization rate 
used here is 3xl0~ 17 s , an appropriate value for high-mass 
star forming regions (e.g. van der Tak & van Dishoeck 2000). 

The model has been run for different temperature values 
(T^gas = ?dust from 10 to 100 K) and volume density fixed 
at the average value of the objects in our sample (n(H2) = 
6 x 10 5 cirT 3 ) . All models start with undepleted abundances 
of the neutral species and run for a time given by the free-fall 
time appropriate for the chosen density (% = 5xl0 4 yr), which 
implies that chemical and dynamical timescales are assumed to 
be comparable. 

Figure[3]shows the model predictions for Df rac , i.e. the deu- 
terium fractionation in species such as N2H + , and the CO de- 
pletion factor /d versus the gas (dust) temperature. The thin 
curves show models where the Gerlich et al. (2002) rate co- 
efficients for the proton-deuteron exchange reactions are used, 
whereas the dotted curve is for models with the (factor of 3) 
larger rate coefficients typically used in chemical models (e.g. 
Roberts et al. 2003, 2004; but see Walmsley et al. 2004; Flower 
et al. 2005 1. In fact, the recently measured rate coefficients still 
suffer some uncertainties, especially due to the importance of 
the so-called "back reactions" of the deuterated forms of HJ 



with ortho-H2, which lower the deuteration depending on the 
unknown ortho/para H2 ratio. 

The data points in Figure [5] suggest that the model is in 
quite good agreement with observations for the majority of the 
objects, despite of the simplicity of it. Considering the uncer- 
tainties in the rate coefficients and in the gas/dust temperature, 
the observed deuterium fractionation is well reproduced in the 
observed range of temperatures. In the case of the CO depletion 
factor, there are three objects that show too large /d for their 
adopted temperatures: IRAS 20126+4107, IRAS 18144-1723, 
and IRAS 19092+0841. In the case of IRAS 20126+4107 and 
IRAS 19092+0841 we note that the excitation temperature de- 
rived from N2H + (1-0) is indeed close to 10 K, suggesting that 
in these cases the deuterated gas may be confined in smaller and 
colder clumps, compared to the high mass star forming region 
where they are embedded. Therefore, lower dust/gas tempera- 
tures (close to 10 K) should be more appropriate. 

However, for IRAS 18144-1723, the N 2 H + (1-0) excitation 
temperature is quite large ( 26 K), so that some other mecha- 
nisms may be present to maintain a large fraction of CO frozen 
onto dust grain (unless the dust temperature is at least 5 K lower 
than the gas temperature). One possibility could be the trapping 
of CO molecules in H2O ice. Indeed, if a large fraction of CO 
is trapped in water ice, dust temperatures between 30 and 70 K 
are needed in order to release all the CO back in the gas phase 
(Collings et al. l2003l Viti et al.|2004i. To test this possibility, 
we run models where a certain fraction of CO molecules were 
assumed to stick onto dust grains with binding energies appro- 
priate for H20 ice (4820 K; Sandford & Allamandola 1990). 
In the case of IRAS 18144-1723, to reproduce the observed 
/d value (3.8) at the measured temperature (23.6 K), 70% of 
solid CO needs to be trapped into water and this is what is 
shown in Figure [3]by the dashed curves, which is also able to 
reproduce the / D value observed in the warmest source (IRAS 
19092+0841). However, the CO trapping does not affect the 
deuterium fractionation, which is limited by the gas tempera- 
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ture. We then conclude that a significant fraction of CO may be 
indeed trapped in H2O ice, although more accurate determina- 
tion of the dust and gas temperature are needed to give more 
quantitative estimates. 

4.3. Where are the pre-stellar cores? 

The most important result of this work is the detection of NaD + 
emission in 7 of our sources, with values of deuterium fraction- 
ation close to those found in low-mass starless cores. In this 
section we discuss this result, focusing the attention on the ori- 
gin and the nature of the cold and dense gas that gives rise to 
this emission. 

While Crapsi et al. (2005 1 have established that in their ob- 
jects the N2D + emission arises from the densest core nucleus, 
the angular resolution of our data does not allow one to deter- 
mine the accurate location of this emission. In fact, as already 
noted in Sects. l4~Tl and l4. 21 it is well known that high-mass stars 
typically form in clusters, so that the surroundings of a mas- 
sive protostar often show complex morphology (see Molinari 
et al. 2002 , Cesaroni et al. 2003 Fontani et al. 2004a 1 and may 
be fragmented in objects with different masses and in differ- 
ent evolutionary stages (see e.g. Kurtz et al. 2000 Fontani et 
al. 2004b), whose angular separation can be comparable to or 
even smaller than the angular resolution of the observations. 

For this reason, in principle the observed cold and dense 
gas responsible for the N2D + emission may be associated with 
the high-mass protostellar object or with another object located 
very close to it. In the first case, such gas is located in the most 
external shell of the parental cloud not yet heated up by the 
high-mass (proto-)star: a 'record' of the early cold phase. As a 
matter of fact, relatively high values of deuterium fractionation 
have been found also in some hot cores (e.g. Oloffson 1984 
Turner 1990), i.e. the molecular environment in which mas- 
sive stars have been recently formed. Such a scenario would 
indicate that the molecular cloud in which low-mass and high- 
mass stars are born have similar chemical and physical con- 
ditions. In the second case, the N2D + emission is due to one 
or more molecular cores in the pre-stellar phase located close 
to the central protostar. In fact, high-angular resolution obser- 
vations of one of our targets, IRAS 05345+3157, have clearly 
shown that the molecular gas is fragmented into several cores 
(see Molinari et al. 2002), separated on average by ~ 5"- 10". 
Such a separation is smaller than the angular resolution of our 
data. Therefore, the observed N2D + emission might arise from 
some of these cores. Given the comparable distances, in princi- 
ple this could be the case also for the other sources. Of course, 
to solve this problem observations with higher angular resolu- 
tions are absolutely required. 

5. Conclusions 

We have observed several rotational lines of N2H + , N2D + , 
C I7 0, DCO + and sub-mm continuum emission with the IRAM- 
30m telescope and the JCMT in 10 high-mass protostellar can- 
didates. Our main goal is to measure the deuterium fractiona- 
tion through the ratio A^(N 2 D + )/A^(N 2 H + ), and the CO deple- 
tion factor (ratio between 'expected' and observed CO abun- 



dance), in order to shed light on the chemical and physical 
properties of the molecular clouds in which high-mass stars are 
born. The main results of this study are the following: 

- We have detected N2D + emission in 7 out of 10 sources 
of our sample, with values of A^(N 2 D + )/A^(N 2 H + ) between 
~ 0.004 and ~ 0.02, with an average value of 0.015. This 
value is 3 orders of magnitude larger than the interstellar 
value for the relative D/H abundance (~ 10 ), and close to 
that found by Crapsi et al. J2005> in their sample of low- 
mass starless cores. 

- For almost all the observed sources the sub-mm continuum 
map and the maps of the N2H + (1-0) and (3-2) integrated 
line intensity show a similar distribution of the emission, re- 
inforcing the idea that N2H + and sub-mm continuum trace 
the same material. 

- The CO 'expected' abundances are comparable to or 
marginally higher than the values measured from the obser- 
vations in 8 sources, and much higher in 2 sources (IRAS 
20126+4104 and IRAS 20343+4129). The median value 
is 3.2, which is smaller than that found in low-mass pre- 
stellar cores. However, the non-depleted gas along the line- 
of-sight may strongly affect these results, since the source 
distances are larger than 1 kpc. We thus believe that the 
'real' CO depletion factors estimated by us are lower lim- 
its. 

- A simple chemical model suggests that the majority of 
the observed sources have conditions similar to low-mass 
pre-stellar cores (large CO depletion factors and deuterium 
fractionation). The observed deuterium fractionation and 
CO depletion can be reproduced if the gas and dust tem- 
perature are ~ 20 K, i.e. close to that observed by NH3 
or CH3C2H observations and to the N2H + (1-0) excita- 
tion temperature. Although some of the studied sources 
(IRAS 20126+4104 and IRAS 19092+0841) may indeed 
host dense and cold clumps with physical and chemical 
characteristics identical to low-mass pre-stellar cores, on 
average the present observations are consistent with enve- 
lope material recently heated by the central massive star, 
suggesting that the studied sources may have experienced a 
pre-stellar phase with gas and dust temperatures quite sim- 
ilar to those observed in the nearby low-mass cores (~ 10 
K), despite the larger volume densities. It appears that large 
fractions of CO may be trapped in H2O ice, although high 
angular resolution observations are needed to confirm this. 

Our results allow us to conclude that most of our sources 
host cold and dense gas with properties similar to those found 
in low-mass starless cores on the verge of the gravitational col- 
lapse. The problem to solve now is the location and the origin 
of this gas. We propose two scenarios: in the first one, the cold 
gas is distributed in an external shell not yet heated up by the 
high-mass protostellar object, a remnant of the parental mas- 
sive starless core. In the second one, the cold gas is located in 
cold and dense cores close to the high-mass protostar but not 
associated with it. This scenario is supported by high-angular 
resolution observations of some of our objects, which have re- 
vealed their clumpy structure. Clearly, observations at higher 
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Fig. 3. Top left. Deuterium fractionation versus kinetic temperature derived from NH3 observations. The arrows indicate the 
deuterium fractionation upper limits. The curves represent the prediction of the theoretical model described in Sect. 4.2. Solid 
curve is model prediction when using the rate coefficients measured by Gerlich et al. (2002) for the proton-deuteron exchange 
reactions. Dotted curve is from the standard model, with (factor of about 3) larger rate coefficients. Dashed curve is from the 
standard model, when 70% of CO molecules are trapped in H 2 ice. 
Bottom. Same as Top for the integrated CO depletion factor. 

Top right. Deuterium fractionation versus depletion factor for the same models and data in the previous two panels. Note that 
the presence of trapped CO does not significantly affect the deuterium fractionation. 
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Fig.A-1. Spectra of the N 2 H + (1-0) line (left panels) and the N 2 D + (2-1) line (right panels) for IRAS 05345+3157, IRAS 
05373+2349, IRAS 18144+3157, IRAS 1851 1+0146 and IRAS 18517+0437 (from top to bottom). All spectra have been taken 
at the position of the N 2 H + (1-0) line emission peak. 
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Fig.A-5. Same as Fig. lA-3l for IRAS 18144+3157. We do not show the C 17 map because it is too noisy, nor the DCO + map 
because the source is not detected in this tracer. 



Fontani et al.: Massive pre-stellar cores 



21 



1 °51 '00' 



1°50'40' 



1 °50'20' 



1°50'00' 



o 
o 
o 

CN 



o 

CD 
Q 



1°50'40' 



1 o 50'20" - 



+ (i-o)- 




N 2 H + (3-2) 




D 50'40' 



3 50'20 



1 o 50'40' 



1 °50'20' 




1 8 h 53 m 38 s 



1 8 h 53 m 38 s 



R.A.fJ2000 



RA(J2000 



Fig.A-6. Same as Fig.lA^3lfor IRAS 18511+0146. 
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RA(J2000) 

Fig.A-7. Left panel: integrated N 2 H + (1-0) (top panel) and (3-2) (bottom panel) line emission in IRAS 18517+0437. Right 
panel: integrated emission of the C n O (2-1) line. In each panel, the filled triangles at map center indicate the position of the 
1.2 mm emission peak (Beuther et al. 2 002b . 
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Fig.A-8. Same as Fig.lA^3lfor IRAS 19092+0841. 
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Fig.A-9. Same as Fig. IA-7l for IRAS 20126+4104. The triangle corresponds to the position of the 850 fim emission peak 
(Williams et al. l2004l 
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Fig. A-10. Same as Fig.|7v9lfor IRAS 20216+4107. 
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Fig. A-ll. Same as Fig.lA^Jfor IRAS 20343+4129. 
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Fig. A-12. Same as Fig.|A21for IRAS 22172+5549. 



